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This  study  of  propellant  fracture  has  been  subdivided  into  three 
parts:  (l)  pertinent  elastic  fracture  developments,  (2)  viscoelastic 
fracture,  and  (3)  dewetting. 

In  (l)  a  fracture  criterion  is  developed  for  stress  fields  producing 
mixed  mode  crack  propagation  by  drawing  a  comparison  with  fracture  in 
uncracked  geometries.  Such  mixed  fields  are  prevalent  for  cracks  in 
rocket  motors  with  star  shaped  cavities. 

In  (2)  a  means  for  predicting  crack  growth  and  velocity  has  been 
developed  for  linear  viscoelastic  meterials  strictly  on  the  basis  of 
stress  intensity  factors.  A  fracture  characterization  of  critical  stress 
intensity  factor  for  a  CTEB  propellant  has  been  completed,  which  facilitates 
the  application  of  the  theory. 

In  (3)  extensive  dewetting  experiments  have  been  completed  which 
have  confirmed  prior  theoretical  developments  of  a  quasielastic  dewetting 
theory.  Not  only  has  the  concept  of  a  stress  criterion  for  dewetting 
been  substantiated,  but  the  form  of  the  surface  predicted  by  the  theory 
has  been  verified.  y 


FOREWORD 

This  report  summarizes  the  results  of  work  on  a  research  project 
entitled  Propellant  Fracture  Studies,  which  is  funded  by  the  Naval 
Weapons  Center,  China  Lake,  California. 

The  report  is  prepared  in  three  sections.  The  first  section 
contains  ideas  and  observations  concerning  elastic  fracture  in  combined 
symmetric  and  skew  symmetric  stress  fields.  Validation  of  the  theory 
is  made  with  gl&ssy  polymers,  and  the  work  was  performed  as  a  precursor 
to  formulating  similar  results  for  viscoelastic  materials. 

Section  2  deals  directly  with  viscoelastic  crack  propagation  in 
symmetric  modes  and  develops  a  much  more  simplified  theory  than  exists 
in  the  literature  for  handling  viscoelastic  fracture  .  The  author  would 
like  to  acknowledge  many  discussions  of  this  material  with  LCDR  C. 
Hertzler. 

Section  3  deals  with  dewetting,  which  is  the  phenomenon  that  makes 
propellant  different  from  ordinary  viscoelastic  materials.  Extensive 
studies  of  dewetting  have  been  reported  previously  on  this  contract 
and  this  section  continues  the  report  of  those  studies.   Acknowledgement 
is  given  to  LCDR  J.  E  Wood,  who  has  contributed  some  of  the  results  of 
his  PhD  thesis  work,  which  is  on  this  topic,  for  this  report. 
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English  Symbols 


6   6   6 

a  ,b  ,c   =  elastic  boundary  quantities 


V   V   V 

a  ,b  3c  =  viscoelastic  boundary  quantities 

b  =  strip  width 

B,B- ,Bp  =  definition  of  boundary 

c  =  crack  length 

E  =  modulus 

f,g,F,G,$  =  general  functions 

G  =  shear  modulus 

h  =  heaviside  step  function 

I, ,Ip,Jp  =  stress  invariants 

K  =  bulk  modulus 

Kn ,  Kp  =  stress  intensity  factor 

r  =  polar  coordinates 

t..  =  time  of  arrival  of  crack 

t  =  sheet  thickness  or  time 

U .  =  displacement 

u  =  strain  energy 

x,y  =  cartesian  coordinates 


Greek  Symbols 


=  experimental  constant 

=  crack  orientation  angle 

g 

ij      =  strain 

^        =  Geometric  factor 

v        =  Poissons  ratio 

$       =  stress  function 


a.  .  =  stress 

a,cr  ,ff   =  applied  stress  at  the  boundary 
n  s 


9       =  polar  coordinate 
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SECTION  1 

SOME  OBSERVATIONS  ON  ELASTIC  FRACTURE  UNDER 
COMBINED  LOADING 


INTRODUCTION 

Much  of  the  traditional  theory  of  fracture  has  been  cast  in  terms 
of  failure  surfaces  in  principal  stress  space,  which  have  been  given 
physically  descriptive  names  such  as  maximum  principal  stress,  distortional 
energy,  etc,  and  these  theories  have  been  applied  to  bodies  that  are 
assumed  to  contain  no  dominant  flaws.  On  the  other  hand,  in  more  recent 
times  an  equally  strong  theory  in  the  form  of  "fracture  mechanics"  has 
been  developed  upon  the  premise  that  a  crack  already  exists  within  the 
body.  As  a  result  of  these  two  developments,  fracture  analysis  is 
compartmentalized  into  problems  of  one  type  or  the  other. 

The  two  approaches  have  been  further  dissimiliar  in  the  types  of 
stress  fields  typically  treated.  While  the  use  of  combined  stress 
states  is  very  common  in  the  instance  of  failure  surface  theories,  it 
is  not  at  all  common  to  apply  fracture  mechanics  to  combined  loading 
conditions,  and  in  fact  a  general  criterion  is  not  available  where  more 
than  one  mode  of  fracture  is  present.  In  this  report  an  attempt  is  made 
to  bridge  between  the  two  for  elastic  materials,  linking  them  more 
closely  together  and  providing  a  more  unified  perspective.  By  bringing 
both  of  these  theories  together,  a  general  fracture  criterion  has 
emerged  for  combined  loading  and  mixed  mode  crack  propagation. 


Work  on  the  problem  of  crack  propagation  in  sheets  under  combined  loads 
was  undertaken  by  Erdogan  and  Sih   ,  who  proposed  an  explanation  of  crack 
trajectories  in  plane  two  dimensional  stress  fields  based  upon  experimental 
results  obtained  with  FMM  (polymethylmethacrylate).   Later  Panasyuk,  et  al 

independently  proposed  similar  explanations  based  upon  their  work  with  glass. 

[3] 
Cottrell   ,  has  examined  the  energy  released  from  a  system  when  a  crack 

extends  radially  along  an  arbitrary  trajectory.  He  found  that  the  direction 

of  propagation  hypothesized  by  the  above  investigators  corresponds  to  the 

direction  of  maximum  energy  release,  which  is  consistent  with  the  general 

fracture  theory  proposed  by  Griffith. 

In  the  experimental  studies  that  were  conducted  relative  to  the  before 

mentioned  investigations,  the  stress  conditions  which  precipitated  crack 

instability  were  only  of  secondary  interest;  consequently  observations  are 

offered  herein  to  complete  the  picture  by  describing  the  local  conditions  of 

fracture  at  a  crack  tip  in  an  arbitrary  two  dimensional  stress  field  by  relating 

them  to  fracture  conditions  in  unflawed  geometries.  This  leads  quite  naturally 

to  a  fracture  criterion  representable  in  K,-!^  space.  From  this  development 

further  observations  can  be  made  regarding  the  changing  nature  of  the  stress 

field  in  the  vicinity  of  a  crack  as  it  continues  to  propagate  under  combined 

loading.  Verification  of  the  theory  has  been  demonstrated  on  FMM  for  two 

dimensional  plane  geometries. 

Interest  in  the  combined  loading  fracture  conditions  arises  for  solid  propellam 

rocket  grains  since  many  of  the  cracks  in  rocket  motors  are  subjected  to  both 

symmetric  and  anti  symmetric  loads. 


CORRELATION  OF  FRACTURE  CRITERIA 

For  the  unflawed  case,  i.e.,  for  specimens  without  cracks,  fracture 
in  multiaxial  stress  states  is  describable  in  terms  of  a  surface  in 
principal  stress  space.  Over  the  years,  it  has  been  found  by  experience 
that  the  continuous  surfaces  are  quadric  surfaces  and  are  representable 
by  second  order  equations  in  stress.  A  very  general  representation  of 
such  surfaces  can  be  made  in  terms  of  the  first  two  stress  invariants 
I..  and  Ip;  however  it  is  physically  more  meaningful  to  use  I.,  and  Jp, 
the  second  deviatoric  stress  invariant.  The  equation  of  the  surface 
in  functional  form  is 

f  (1^,  j2)  =  0  (1.1) 

2 
where  I   is  selected  for  dimensional  homogeniety. 

Mathematical  expressions  such  as  equation  (l«l)<  represent  the  manner 
in  which  a  material  reacts  fracturewise  to  applied  loads  and  theoretically  will 
predict  fracture  regardless  of  the  geometry  or  character  of  the  stress 
state.  A  crack  in  a  body  introduces  a  geometrical  change,  but  it 
should  not  influence  the  fundamental  fracture  behavior.  Therefore  a 
similar  function  should  exist,  which  is  related  to  this  one,  and  which 
describes  fracture  of  cracked  bodies. 

For  geometries  containing  cracks,  Barenblatt    has  proposed  that 
the  local  distribution  of  cohesive  forces  over  the  crack  surface  at 
fracture  is  a  material  property,  independent  of  external  loading 
and  body  geometry;  thus  for  any  body  constructed  of  a  given  material 
there  exists  a  function  of  stress  intensity  factors  which  describes 


elastic  fracture  conditions.  For  two  dimensional  stress  states  limited 
to  the  opening  and  sliding  mode,  his  criterion  reduces  to 

*  ^a'  "K2a)  =  °  (1'2) 

where  K_  and  K_  are  the  portion  of  the  stress  intensity  factors  which 
are  due  to  the  cohesive  forces  at  the  crack  tip.  The  total  stress 
intensity  factor  is  obtained  by  adding  the  cohesive  force  contribution 
to  K,  and  IC,  which  result  from  the  external  forces. 

M 

By  Barenblatt ' s L  J  hypothesis  of  no  singularity  at  the  crack  tip 

Thus  from  equations (1.2)  and  (l.U) 

$  (K^,  Kg)  =  0  (1.5) 

To  relate  f  in  equation  (l.l)  to  $  in  equation  (1-5),  §(K ,,  K  )   will 
be  expressed  in  terms  of  invariants.  The  stress  field  in  the  vicinity 
of  a  crackjwhen  cohesive  forces  are  not  considered  ,may  be  written  as 

^  aij  =  ^L  Fij  (9)  +  *2   Gij  (9)      i,J  =1'2       (1'6) 

From  this  equation,  one  can  solve  for  K_  and  K_  in  terms  of  the  two 
in-plane  principle  stresses  or  in  terms  of  the  two  stress  invariants, 
both  of  which  can  be  used  to  relate  fracture  of  cracked  geometries  to 
fracture  of  uncracked  geometries .  Continuing  in  terms  of  invariants , 
equation  (1L.6)  can  be  used  to  formulate  I,  and  J?  . 


r  Ij2  =fx  Q^,  Kg,  9)  (1.7a) 

.  r  J2  =f2  (1^,  I^,  0)  (1.7b) 

Solving  for  IC  and.  K-  and  inserting  them  into  (1.5)  produces  a  description 
of  the  local  conditions  governing  fracture  at  a  pre-existant  crack 

$  (rJ2,  rl^)  =  0  (1.8) 

This  describes  the  conditions  of  crack  propagation  in  some  direction, 

9  >of  the  maximum  tensile  stress,  as  referenced  from  the  original  crack 
o 

orientation.  It  disregards  the  cohesive  forces  and  expresses  the 
fracture  condition  in  terms  of  the  externally  applied  loads  only  just 
as  the  fracture  surface  does.  On  the  other  hand,  the  success  of  the 
fracture  surface  concept  over  a  wide  variety  of  stress  states  and  "body 
shapes  makes  evident  the  fact  that  the  fracture  behavior  of  an  elastic 
material  is  independent  of  the  geometry  of  the  body  of  interest,  and 
it  will  not  be  altered  in  a  fundamental  way  by  the  introduction  of 
different  geometric  factors,  such  as  cracks.  It  seems,  therefore, 
appropriate  to  offer  as  a  hypothesis  that,  when  only  external  forces 
are  considered,  "Fracture  will  occur  in  cracked  geometries  at  the  same 
local  stress  conditions  which  produce  fracture  in  uncracked  geometries." 
In  other  words,  the  crack  will  extend  when  the  invariants  of  stress  on 
the  trajectory  satisfy  equation (lJ.)  and  equation  (1.8).  This  can  be 
expressed  mathematically  as 

*  =  f  (1.9) 


In  other  words,  if  the  fracture  surface  f  is  known,  the  fracture 

condition,  §,  for  bodies  containing  cracks  may  he  obtained  by  replacing 

2  2 

J2  and  I,  with  r  Jp  and  rL.  .  From  this  a  general  fracture  criterion 

can  be  generated  in  terms  of  K,  and  Kp. 


EXAMPLE  PLEXIGLASS 

Figure  1-1  shows  the  uncracked  multiaxial  fracture  behavior  of 

T6 1 
PMM  obtained  by  Thorkildsen  and  Olszewski   ,  where  it  has  been 

plotted  in  two  dimensional  principal  stress  space.  A  more  precise 

reduction  of  this  data  in  terms  of  invariants  makes  possible  the 

correlation  with  results  obtained  by  Erdogan  and  Sih    on  cracked 

specimens  made  from  the  same  material. 

Stress  invariants,  I,  and  I?,  are  defined  as: 

L,  *  a,  +  a2  +  o-s  (1.10a) 

I2  =  -(a1  a2  +  a2  a~  +  o^  cO  (1.10b) 

2 
A  plot  was  made  in  Figure >2  of  Ip  vs  I   for  the  experimental  data, 

which  appears  to  yield  a  clearly  defined  linear  relationship.  The 

data  was  fitted  by  least  squares  to  a  linear  equation,  which  resulted 

in 

.Ix2  =  -2.93I2  +  0.96au2  (1.11) 

where  a  is  the  failure  stress  in  uniaxial  tension.  This  expression 
u 

in  terms  of  invariants  is  postulated  to  be  the  description  of  a 
failure  surface  for  all  possible  tensile  stress  states.   Instead 
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FIG.   1-1    FRACTURE        STRESS      FOR      POLY   (methylmethacrylate) 
(Thorkildsen      and       Olszewski*6*  ) 
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FIG.  1-2     LEAST      SQUARES      FIT     TO     FRACTURE 
DATA     FOR      PMM 
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of  leaving  the  final  expression  in  terms  of  I  ,  which  has  no  particular 
physical  interpretation,  it  is  preferred  to  use  J_,  which  is  related 
to  I1  and  Ip  by 

J2  =  —  +  h  d-12) 

Using  (1JK)  and  (1.12)  an  expression  for  the  failure  surface  of  PMM  is 
found  to  be  of  the  form, 

f(l12,  J2)  =  J2  +  0.0071^  -  Const  =  0  (1.13) 

As  can  be  seen,  there  is  a  very  slight  dependency  on  hydrostatic 
pressure,  but  shear  is  the  predominant  factor  for  fracture  in  tension 
of  PMM. 

As  stated  previously,  it  is  proposed  that  this  criterion  also 
governs  the  onset  of  crack  propagation  for  bodies  containing  residual 
cracks  in  the  same  type  of  material .  In  other  words ,  for  bodies 
containing  a  crack,  fracture  will  occur  by  an  incremental  extension 
of  the  crack  along  the  path  normal  to  the  greatest  tension  when  the 
stresses  on  that  trajectory  combine  to  satisfy  the  fracture  criterion 
found  for  uncracked  specimens .  This  requires  that  the  direction  of 
propagation,  S  ,  be  known  a  priori;  however  this  information  can 
be  calculated  from  the  work  of  Erdogan  and  Sih.  By  equations  (1.8) y    (1.9) 
and  (1.13),  propagation  of  a  crack  is  given  by 

rJ2  +  .007  rl,   -  Const  =  0  (l.l^) 


To  apply  this  criterion,  the  crack  tip  stress  field  is  used. 

The  coordinate  description  is  standard,  and  the  equations  are  classical 

[71 

as  given  by  Sih  and  Liebowitz   . 

a  =  ^-r  Cos  |  [(3  -  Cos  9)K_  +  (3  Cos  9  -  1)1^  Tan  f]       (1.15a) 

r       2(2r)2  2  1  -2  2 

3  =  ^-i  Cos  |  [(1  +  Cos   9)1^  -  3^  Sin  9]  (1.15b) 


o«  = 


rQ=      ;  \i  Cos  |  [lysine  +  KgO  Cos    6  -  l)]  (1.15c) 


T  -= 


2(2r) 


J  in  polar  coordinates  is  given  by, 

J0  =  i  [a2  -  a   aa  +  a2  +   3t2  ]  (1.16) 

2   3   r    r  9    9    r9 

which  in  terms  of  stress   intensity  factors  from  equations  (1.15)  give 

r J2  =  \     {I  ^  +  k  COS0  "  3  C°S  29^  4  +  ^2  Sin  29  "  Sin  9^KiK2 

+  |  [7  -  Cos  9  -  |  Sin2   9]K^}  (1.17) 

From  equations  (115),    I,    can  also  be  readily  computed, 

r  I2  =  2   [K^  Cos2  |  -  JLfe  Sin9  +  K^Sin2  |]  (1.18) 

Combining  equations (1.17)  and  (1J^      for  a  surface  of  the  general  form 
that  was  exhibited  by  PMM  in  equation  (1-1%  we  have 

r(J2  +  ol2)  =  g-  {[|  (7+*+  Cos9  -  3  Cos  29)  +  6a  (l  +  Cos   9)]  K^ 

,+  [    (|  Sin  29  -  Sin9)   -  12  a  Sin9]  K^  (1.19) 

+  [|  (7  -  Cos9  -  |  Sin29)  +  6a(l  -  Cos9)]  K?J 
10 


The  constant  in  (1.1*+)  is  found  by  examining  the  fracture  condition 
of  symmetric  loading,  for  which  IC,  -»  }L    3   IC  =  0  and  9  =0.  From  (1.19) 

[r(J2  +  al2)]  =  Const  =  |(l+12a)  K^2  (1.20) 

Thus  the  criterion  of  fracture  for  arbitrary  combined  loads  is  given 
by 

(1  +  12a)K^c  =  [(^  +  6a)  +  (|  +  6a)  Cos  6  -  |  Cos  2 el  K^ 

+  [|  Sin  20  -  (l+12a)  Sine]^^  (1.21) 

+  [(I  +  6a)  -  (|  +  6a)  CosG  -  £  Sin2e]  K^ 

where  for  PMM,  a  =  0.007- 

To  limit  check  the  skew- symmetric  fracture  condition,  K_  =  0, 
Kp  -♦  Kp  and  Erdogan  and  Sih's  theory  and  measurements  give  6  =  -70.5  . 
Equation  0-21)  gives  the  critical  K?  value  in  terms  of  the  critical 
K_  value, 

_  r  3  a  +  1200-1*  (1 22) 

*2C  ~  L  h   (i  +  3a)  J   ^c  U*22j 

Substantiation  of  this  prediction  from  the  theory  is  provided  by  the 
experimental  results  of  Erdogan  and  Sih,  who  found  for  PMM,  YL-,   =  ^72 
lb  in"3/  and  K^  =  U22.1  lb  in"3/2.  Using  this  value  of  K„  in  (1.22), 
and  a  =  .007  from  ( 1.1^)  the  theory  predicts  IC  =  1+21.5  lb  in"3'  . 
This  very  close  correlation  provides  meaningful  substantiation  to  the 

» 

theory. 


11 


Further  evaluation  of  the  theory  can  be  obtained  by  plotting 
equation  (1.3.)  for  a  variety  of  combined  loading  conditions .  For  any 
given  angle  of  propagation,  9  ,  the  fracture  equation  (1.3)  is  an 
inclined  ellipse  in  L,  K_  space.  However  in  YL  ,  IC  space,  each 
set  of  coordinates  dictates  a  different  loading  condition,  which 
requires  the  angle  of  propagation  to  be  different,  and  the  coefficients 
in  the  fracture  equation  vary,  deviating  from  an  ellipse.  A  plot  of 
the  equation  is  given  in  figure  1-3,  where  it  is  compared  with  the  data 
of  Erdogan  and  Sih  for  a  variety  of  combined  loading  conditions 
obtained  by  using  the  geometry  of  a  crack  inclined  to  the  direction 
of  applied  loads  for  several  different  crack  orientation  angles,  (3. 

As  can  be  seen,  the  theory  fits  the  data  quite  well  for  strongly 
skew-symmetric  loadings,  but  deviates  somewhat  in  some  of  the  loadings 
that  are  predominantly  symmetric .  Of  course  for  those  loadings  that 
are  purely  symmetric,  the  criterion  matches  identically  due  to  the 
fact  that  the  governing  equation  was  matched  to  this  condition. 

SUBSEQUENT  BEHAVIOR 

After  the  crack  has  begun  to  run  under  combined  loads,  the 
geometry  will  appear  as  depicted  in  figure  1^4 for  small  extensions. 
A  solution  of  this  problem  has  been  obtained  by  Andersson   ,  wherein 
he  found  that  the  stress  intensity  factors  K.  and  K~  are 

/\         r  tt  -  9  _  2 

h  *  °  F  L  TTTl        sin  p  Sin  (p  +  eo} 

o 

eo  (1.23) 

%  =  CT  ^  [  I  I  9°]         Sin  P  C°S   (3  +  9o) 

o 
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FIG.  1-3   K,   VERSUS   K2   AT  THE   BEGINNING  OF 
CRACK   EXTENSION   IN   A   CRACKED   PLATE 
UNDER   PLANE   LOADING. 
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FIG.  1-4   GEOMETRY   OF   CRACK   AFTER   EXTENSION 
UNDER   COMBINED  LOADS. 
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However  before  the  crack  extended,  the  stress  intensity  factors  were 
those  of  an  inclined  crack  in  an  infinite  sheet. 

p 
1^  =  cr^/6"  Sin  p       Y^  -   a^/c  Sin  p  Cos  p  (1.2*0 

When  these  expressions  are  substituted  into  (1.23)  it  reduces  to 

% 

h  =  [  ^rr]      [Ki Cos  eo  -  h  Sin  ^ 

o 


3»  -  [  ¥TT  ]       ^  Sin  6o  "  "a  Cos  eo^ 

O 

This  set  of  equations  suggests  that,  in  effect,  propagation  can  be 
considered  as  a  rotation  of  the  K.K  axis  system,  and  from  the 
experimental  evidence  available,  the  rotation  is  such  that  it  reduces 
the  value  of  1L.  In  other  words  the  crack,  operating  under  the 
observed  laws  of  behavior  mentioned  previously,  will  propagate  in 
the  direction  of  the  maximum  tension  and  this  will  reduce  the  value 
of  Kp.  Furthermore  by  extrapolation,  subsequent  increments  of 
propagation  will  be  along  the  direction  of  instantaneous  maximum 
tension,  which  will  continue  to  reduce  Kp,  driving  it  to  zero  so 
that  ultimately  propagation  continues  in  the  direction  normal  to 
the  maximum  tensile  field  stress  of  the  body. 

To  illustrate  this,  consider  figure  1-4 as  an  example, where  1  =  0 
and  p  =  70  .  The  angle  of  maximum  tension  is  9  =  -30  ,  or  -tt/6. 
The  corresponding  point  on  the  fracture  envelope  in  K, ,  Kp  space  from 
(1.21) and (1.24)  is  ^  =  UU2  lb  in"3/2  and  Kg  =  155  lb  in"3/2.  After 
propagating  an  infinitesimal  amount,  X,  the  new  stress  intensity 

Ik 


(1.25) 


factors  from  (1.25)  are  1^  -  UU7.U  lb  in"3/2  and  Kp  -  81*-. 3  lb  in"3/2. 

Thus -in  the  first  infinitesimal  increment  of  propagation,  Kp  is 

reduced  by  1+6$.  This  change  in  K_,  which  produces  a  change  in  the 

loading  that  the  crack  sees,  causes  a  change  in  the  stress  field. 

These  changes  reduce  Kp  further,  and  this  interaction  continues  until 

Kp  is  nonexistent  and  the  crack  is  propagating  in  a  symmetrically 

loaded  condition. 

Another  interesting  way  of  looking  at  equation  (1-25)  is  to  take) 

the  square  root  of  the  sum  of  squares .  This  results  in 

1  1 

/yp    AO   2  p      p   2 

[££  +  to     =  f(e0)[<  +  sjj]  (1.26) 

In  K,j  Kp  space  the  square  root  of  the  sum  of  squares  is  the  length 
of  a  position  vector.  Since  f(9  )  £  1,  this  equation  shows  that  the 
position  vector  to  the  fracture  envelope  in  K  space  is  becoming  shorter, 
or  that  the  "resultant"  stress  intensity  factor  is  growing  smaller  as 
the  crack  propagates.   If  it  decreases  by  a  sufficient  amount,  the 
crack  should  stop.  The  maximum  reduction  for  the  first  increment 
occurs  at  the  largest  0  ,  which  is  produced  by  pure  shear  loading. 
In  that  instance  0  =  -70.5°,  for  which  f(0  )  =  0.85 .  This  turns 
out  to  be  enough  to  stop  propagation  of  the  crack  after  it  has  begun. 
In  the  experimental  results  of  Erdogan  and  Sih  on  PMM  in  pure  shear, 
the  crack  subjected  to  pure  shear  stopped  after  propagating  a  short 
distance;  however,  this  was  the  only  crack  geometry  studied  for  which 
the  crack  stopped  close  to  where  it  had  begun  to  propagate. 
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These  observations  have  been  discussed  within  the  context  of 
experimental  data  on  PMM  and  glass.  Similar  data  on  inclined  cracks 
in  plates  are  not  available  for  metals.  If  they  were,  the  validity 
of  the  observations  could  be  checked  further  and  the  influence  of 
plastic  yielding  could  also  be  studied;  however  preparation  of  metallic 
sheet  specimens  with  real  cracks  of  arbitrary  angles  of  inclination 
is  quite  difficult.  Meanwhile  the  available  evidence  on  two  elastic 
materials  appears  to  substantiate  a  correlation  between  the  two 
failure  domains . 

SUMMARY 

Observations  have  been  made  which  link  fracture  criteria  in 
unflawed  specimens  to  crack  instability  in  bodies  containing  residual 
cracks s  and  from  this  a  criterion  for  fracture  under  arbitrary  loading 
is  postulated  in  terms  of  an  envelope  in  stress  intensity  factor  space. 
Furthermore  an  explanation  is  offered  for  crack  behavior  subsequent 
to  the  onset  of  propagation  which  is  based  upon  the  concepts  generated 
concerning  the  surface  in  stress  intensity  factor  space.  Experimental 
results  for  two  materials,  glass  and  polymethylmethoerylate ,  verify 
the  theory  for  two  dimensional,  combined  loadings  involving  the 
opening  mode  and  the  sliding  mode  of  crack  propagation. 
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SECTION  2 
VISCOELASTIC  FRACTURE 


SECTION  2 
VISCOELASTIC  FRACTURE 
Introduction 

In  the  last  three  years  viscoelastic  fracture  theory  has  been 
advanced  along  much  more  fruitful  avenues  than  those  which  the  previous 
efforts  produced.  By  applying  principles  of  conservation  of  energy  at  the 
crack  tip,  Knauss   has  shown  that  a  differential  equation  governing 
crack  length  history  can  be  generated.  Several  of  his  works  continue  this 
development, and  in  each  case  an  unknown  length  parameter  manifests  itself 
in  the  differential  equation.  The  parameter,  which  must  be  determined 
experimentally,  seems  to  be  inherent  to  the  method  of  solution. 

Since  this  parameter  is  a  bother  to  deal  with  and  since  it  is 
artificial  in  the  sense  that  there  is  no  physical  basis  for  its  existence, 
it  was  felt  that  a  worthwhile  contribution  could  be  made  by  devising  a 
means  of  formulating  viscoelastic  fracture  without  the  parameter.  The 
means  by  which  this  was  finally  done  was  through  the  viscoelastic  correspond- 
ence principle. 

The  classical  correspondence  principle  normally  applied  in  the 
solution  of  boundary  value  problems  in  the  linear,  quasi-static  theory  of 
viscoelastic,  does  not  apply  for  crack  propagation  problems.  In  that 
instance  new  boundary  is  being  created  as  the  crack  propagates,  and  the 
conditions  on  the  boundary  are  changing  in  type;  i.e.  they  are  changing 
from  displacement  to  stress  boundary  conditions  and  vice  versa.  In  this 
kind  of  problem  there  will  be  points  of  the  boundary  where  only  partial 
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histories  of  the  boundary  conditions  are  known;  consequently  the  integral 
Laplace  Transform  cannot  he  computed,  and  the  correspondence  principle  is 

not  applicable. 

[o] 
Graham    has  extended  the  correspondence  principle  to  encompass 

problems  of  this  type.   (An  explanation  of  his  development  is  given  in 

Appendix  A.)  By  using  Graham's  method  of  solution,  two  very  large 

categories  of  problems  have  been  solved.  Linking  these  solutions  with 

Knauss'  experimental  findings  and  theoretical  concepts  has  lead  to  a 

method  of  viscoelastic  fracture  analysis  based  upon  stress  intensity  factors, 

The  method  parallels  elastic  fracture  theory  very  closely  and  combines  the 

work  of  other  investigators  in  such  a  way  that  the  length  parameter  is 

no  longer  needed. 

Stress  Analysis 

One  very  common  class  of  crack  problems  for  which  the  extended 
correspondence  principle  is  applicable  is  shown  in  Figure  2-1,  where 

the  lower  portion  of  the  body  has  been  removed  and  replaced  with 

T27l 
appropriate  boundary  conditions.  The  two  solutions  by  Graham     in 

the  literature  that  use  the  extended  correspondence  principle  are  for 

bodies  of  infinite  extent.  However  the  method  is  applicable  to  finite 

geometries.   In  the  case  of  stress  boundary  conditions,  like  those 

pictured  in  Figure  2-1,  the  application  is  quite  simple.   In  the  case  of 

displacement  boundary  conditions,  which  is  treated  later,  the  boundary 

conditions  have  to  be  recast  into  an  acceptable  form;  therefore 
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Crack   Surface 

' w«. 


BtCt)- 


Figure  2-1  Definition  of  boundaries  for  the 
Correspondence  Principle 


the  boundaries  may  be  finite  or  infinite.  They  are  loaded  with 

normal  stresses  and  are  free  of  shear  stress.  The  stress  distributions 

are  symmetric  so  that  displacements  on  the  imaginary  extension  in  front 

of  the  crack  are  zero.  The  boundary  conditions  vary  with  time 

according  to  the  following  description: 

a     =  0  on  B 

s 


a  =  B(x,t)      on  Bn(t) 
n    v  >  '         lv  ' 


(2.1) 


u  =  0 
n 


on  B2(t) 


These  conditions  correspond  to  column  2  of  the  matrix  given  in  equation 
(A.l)  of  Appendix  A. 

The  extension  of  the  correspondence  principle  by  Graham  permits  the 
construction  of  the  viscoelastic  solution  for  a  moving  crack  from  an 
elastic  solution  for  a  moving  crack.  This  in  turn  may  be  generated  from 
the  stresses  and  displacements  of  a  stationary  crack.  For  an  incompressible 
elastic  body,  the  field  quantities  in  the  neighborhood  of  the  crack  tip 
are  given  by 

(2.2a) 


a..  =  -i-  f.  .(9) 
1J  4^    1J 
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3Kl  /2r~ 

ui  =  5E-VTrgi(e) 


(2.2b) 

The  crack  is  now  Imagined  to  propagate  slowly.  For  most  elastomers 
and  all  propellants,  the  velocity  of  propagation  is  quite  low,  which 
makes  the  inertial  terms  in  the  equations  of  motion  small  and  neglectable . 
Therefore  the  stress  and  displacement  fields  are  the  same  as  those  given 
in  equation  (2)  except  that  the  coordinate  system  to  which  the  above 
quantities  are  referenced  is  now  moving  since  it  is  attached  to  the 


crack  tip.  The  dynamic,  elastic  solution  is 

K?(t) 

<*-(*)  =  ~ f.  Ie(t)] 

1J  VifTrTtl     1J 


(2.3a) 


.  K®(t)  /TTTT 


(2.3b) 


i*  '   8  "  E   V  tt 

It  is  a  simple  matter  to  reference  these  quantities  to  an  inertial 
coordinate  system.  As  stated  previously  the  loading  and  geometry  is 
taken  such  that  only  the  opening  mode  of  propagation  is  present;  there- 
fore 

r(t)  cos9(t)+c(t)  =  x  (2.4) 


e(t) 


Figure  2-2  Inertial  Coordinate  System 
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Equations  (2.3)  for  the  dynamic  elastic  problem  became 

K*(t) 


*;<(t)  = 


F,,[e(t)l 


V2n(x-c(t))   1J 


1 


(2.Ua) 


(2.Ub) 


To  use  the  correspondence  principle,  the  transform  of  equation  (2.k) 
is  needed. 

* 

e* 


K*(t) 


^   Lv^x-cft)]    1J 


F,,[e(t)] 


By  equations  (k.k)   of  Appendix  A, 

K*(t) 


v* 

J.  . 
1J 


V2n[x-c(t)] 


F..[9(t)] 


E-E* 
v-*  v* 


(2.5a) 


(2.5b) 


(2.6a) 


u. 


e* 


3E* 


Ke(t)  /2£x^t)] 


G.[e(t)] 


E  -»E* 


For  elastic  stress  intensity  factors  in  -which  there  are  no  material 
properties,  equation  (2.6a)  is  readily  inverted  and  yields 

aV.(t)  =a?.  (t) 

From  this  result,  by  comparison  of  the  terms, 

Kj(t)  *,lj<t) 
where  the  stress  intensity  factor  for  viscoelastic  materials  is 
defined  to  be  the  same  as  for  the  elastic  case.   Indeed  when  the 


(2.6b) 


(2.7) 


(2.8) 
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boundary  conditions  are  all  in  terms  of  stress,  the  two  are  identical, 
which  indicates  that  the  definition  chosen  is  consistent.  In  other 
words,  we  are  here  defining  the  viscoelastic  stress  intensity  factor  as 

f  Ce(t)] 

that  factor  which  when  multiplied  by    "         gives  the  stress 

V2n(x-c(t)) 
components  in  the  vicinity  of  the  crack  tip.  This  definition 

possesses  the  appropriate  form  as  one  proceeds  to  the  limit  of  elastic 

materials.  Also  the  stress  intensity  factor,  contains  information 

concerning  geometry  and  boundary  conditions,  as  it  does  in  the  case  of 

elastic  materials. 

For  many  geometries  fitting  the  boundary  conditions  of  equation 

(2.1),  the  elastic  stress  intensity  factor  is 

K-,_  =  \oVttc  (2.9) 

where  \  is  a  geometric  modifiying  factor  that  can  be  independent  of  time, 
but  is  often  a  function  of  c .  If  a  is  the  applied  stress  at  the  boundary, 
then 

K^(t)  =  K^(t)  =  X[c(t)Kt)^FftT  (2.10) 

Fracture  Postulate 

Fracture  studies  in  recent  ye  ars  have  led  to  the  conclusion  that 
flaws  in  materials  will  grow  when  under  load.  This  continues  until  the 
flaw  becomes  of  critical  size, at  which  time  the  specimen  undergoes 
unstable  crack  propagation  and  fracture.  This  has  been  found  to  be 
true  for  metals,  especially  in  fatigue,  and  for  nonmetals  alike. 

Knauss     has  published  theoretical  developments  of  these  ideas 
for  viscoelastic  materials  and  has  shown  that  the  important  parameter 
of  interest  for  rate  sensitive  materials  is  crack  velocity.  In  the 
elastic  case,  the  threshold  of  rapid  crack  propagation  is  often  the 
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focus  of  interest;  however  in  viscoelastic  materials  the  crack  is 
found  to  propagate  for  almost  any  load  level.  The  velocities  are 
infinites imally  small  for  low  loads  and  increase  as  the  body  is  more 
highly  stressed.  The  question  to  be  answered  then  becomes,  not  at  what 
load  level  does  the  crack  propagate  but  what  is  the  crack  velocity  for 
any  given  load.  As  defined,  the  viscoelastic  stress  intensity  factor 
contains  information  relative  to  the  boundary  conditions;  therefore 
•we  desire  the  relationship  between  stress  intensity  factor  and  crack 
velocity. 

This  brings  us  to  the  fracture  postulate: 

L  lx  /JBody  1     lv  /JBody  2 

then  [c(t)]_  ,  .  =  [6(t)]_  ,  Q 
L     Body  1     v   JBody  2 

In  words  this  states  that  for  every  stress  intensity  factor  there  is  a 

corresponding  unique  crack  velocity. 

Knauss     has  performed  experiments  on  Solithane  113  >  a  polyurethane 

[22] 
elastomer,  which  have  substantiated  this  to  a  degree  (Francis  et  al 

have  duplicated  these  tests  and  done  additional  work  on  other  types 

of  loads.)  For  all  temperature  ranges  tested,  the  material  was  elastic 

in  the  sense  that  it  did  not  relax;  however  the  crack  behavior  was 

viscoelastic  and  exhibited  a  velocity  that  varied  with  loading.  Knauss 

postulates  that  in  this  type  of  material  when  the  strain  at  the  crack 

tip  is  the  same  in  two  different  bodies,  the  same  crack  velocity  will 

be  produced.  Of  course  in  a  body  that  is  behaving  elastically  so  far 

as  stress  and  strain  are  concerned,  the  same  strain  level  is  equivalent 

to  the  same  stress  intensity  factor.  Thus  in  the  postulate  of  fracture, 

we  are  following    the  previous  ideas  of  Knauss.  Furthermore,  it 

states  that  if  the  stress  intensity  factor  history  can  be  calculated, 
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the  crack  velocity  equation  can  be  determined  from  the  experimental 
characterization  of  the  material.  Once  this  is  known,  the  crack 
history  can  he  obtained  by  quadrature.   This  will  be  demonstrated  later. 

Fracture  Characterization 

Fracture  characterization  has  been  done  for  only  a  few  materials, 
and  it  has  normally  been  accomplished  with  an  infinite  strip  containing 
a  semi-infinite  crack.   It  could  just  as  well  have  been  done,  however, 
on  the  sheet  with  the  internal  crack  that  was  described  earlier  in  this 
section. 

To  apply  the  correspondence  principle,  we  appeal  to  the  extended 
version  in  Appendix  A.  The  geometry  is  shown  in  Figure  2-3>  where  the 
boundary  is  given  a  normal  displacement  in  the  vertical  direction  on  Bo 


Figure  2-3  Displacement  Boundary  Conditions 


The  actual  boundary  conditions  are  given  below. 


on  B-^t) 

on  B  (t) 

on  B2(t) 

on  B2  (t ) 


uj  =  Un(x,t)  on  B3 
uj=  Ug(x,t)  on  B3 


V 

CTs 

= 

0 

n 

= 

0 

V 

u 
n 

= 

0 

uv 

= 

0 

(2.11) 


2k 


These  must  be  redefined  in  terms  that  will  satisfy  the  conditions  of 

the  extended  correspondence  principle.  This  can  be  accomplished  with 

the  following  set: 

u^  =  Av(x,t)  on  B-^t)  +  B2(t) 

oj  =  0       on  B-^t) 

us  =  0       on  B2(t) 

uj  =  U  (x,t)  on  B~ 

v 
us  =  us(x>"^)  on  B3 

The  normal  displacements  are  specifiable  over  the  crack  surface  without 
material  properties  because  the  far  field  boundary  conditions  are  on 
displacements.  This  now  allows  the  construction  of  an  equivalent 
elastic  problem,  whose  boundary  conditions  are 

u^  =  Ae(x,t)  on  B^t)  +  B2(t) 

a®  =  0  on  B1(t)  +  B2(t) 

u®  =  Un(x,t)  on  B3 

u|  =  Ug(x,t)  on  B3 

Since  the  transforms  are  known  for  all  time, 

u®  =  Ae  (x,p)  on  B^t)  +  B2(t) 
of   =  0        on  B1(t)  +  B2(t) 

y  y 

un  =  Un  (X'P^  on  B3 

y  y 

us  =  Us  ^x,p^  on  B3 
represent  the  boundary  conditions  to  the  corresponding  viscoelastic 
problem.  From  Graham  the  viscoelastic  solution  is  obtained  from 
(A.*0. 
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The  stress  intensity  factors  for  the  strip  are  developed  in  Appendix  B. 
Using  the  results  of  equation  (B.20a)  in  equation  (2.5)  for  v  =  ■§■, 


e* 

r.  . 


E 


«0(t) 


-i  * 


V2n[x-c(t)] 


F. .(0) 


9  =  0 


(2.12a) 


* -<$[•.« 


2Cx-c(t)] 


G.  (tt)6 
i 


x  <  c  ,  6  =  1 

o' 


X 


>  c0,  6=0,(t-tl)   <2-12» 


By  the  correspondence  principle 


v*      ,/lrtj       * 

CTij  =Vr  EREL  P 


e    (t) 

ov    ' 


V2n[x-c(t)]  m 


F.  .(0) 


(2.13a) 


u 


V*    _^3b 


ov    *y  TT 


-c(t?] 


G.(tt)  6 
i 


I  2 

Inverting  the  transforms,  the  viscoelastic  solution  is  obtained. 


(2.13a) 


•0(T) 


V2tt[x-c(t)]  . 


(3T  Fid(0) 


(2.ll»a) 


uT(t)=^eo(t)A/^-£it)l  Gjn)  6 


1"  '      2  O        "\l  TT         1 

It  is  evident  from  these  equations  that  now  the  simpler  expressions 
are  the  displacements,  since  the  boundary  conditions  are  of  that  type. 
The  stresses  are  now  more  complex  due  to  the  material  properties  that 
are  contained  in  the  stress  intensity  factor. 

For  a  step  strain  history,  which  is  most  usually  the  one  that  is 
applied, 

e  (t)  =  e  l(t) 

ov  '         o  v  ' 


(2.lUb) 


(2.15) 
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The  stresses  and  displacements  reduce  to 


o-,  A  (t )  =  \/—  <  -=— =r  +  "e  J  -7-  dT>  F, <  (0 ) 


(2.16a) 


U 


V2rr(x-c  ) 


o  [2nfc-c(T)]] 


3/2 


ij 


^g   ,/2[x-c(t)] 


uv(t)  =  *f-  e  A/*L*-<^/J  G.(tt)6  (2.16b) 

lv  '       2    O  V      TT         1  .' 

The  stress  intensity  factor,  by  comparison  with  equation  (2.Ua),  and 
by  the  previously  given  definition  would  be 

K-.(t)  =A/-=-Cx-c(t)]e  <— —= — -  +  nj   ry-  dT)  (2.17) 

V  3         °(y2n[x-coJ      o  [2tt[x-c(t)]]3/2   J 

This  is  a  very  complicated  expression  ,and  one  not  very  conducive  to  data 
reduction  from  laboratory  results.  However,  a  very  good  approximation 
can  be  made  , which  makes  the  data  reduction  quite  simple. 

A  distance  out  in  front  of  the  crack  in  the  strip  ,the  stress  field 
is  that  of  a  strip  without  a  crack.  If  the  specimen  is  long  enough 
compared  to  the  crack,  a  large  portion  of  the  body  will  be  subject  to 
the  stress  given  in  equation  (B.17).  With  that  approximation,  the  stress 
intensity  factor  is 
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^(t)  =Yf  a(t)  (2.18) 

where  a(t)  is  the  stress  acting  on  the  strip  computed  on  the  basis  of 
the  net  cross -sectional  area,  i.e.  with  the  crack  length  removed.  In 
that  case,  from  equation  (2.5a) 
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e* 

r .  . 
1J 


o-(t) 


V2n[x-c(t)]  _ 


F.  .(0) 


(2.19) 
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Because  there  are  now  no  material  properties  in  the  expression 

v*    e* 

a .  .  =  a . .  and 

oTM   =#      ■     gW  F..(0)  (2.20) 

Therefore  the  viscoelastic  stress   intensity  factor  is 

Ki(t)  =\/¥  a(t)  (2-21) 

Since  the  displacements  are  fixed  at  the  boundary,  the  stresses  in 
the  strip  are  relaxing,  and  K,(t)  will  he  decreasing  in  time  as  the 
crack  propagates.  This  can  be  shown  from  the  viscoelastic  stress - 
strain  relationships  obtainable  from  the  elastic  solution  of  Appendix  B. 
Equation  (B.17)  for  an  incompressible  material  gives 

a  =]f~*   e0  (2.22) 

By  the  regular  correspondence  principle,  the  stress  in  the  viscoelastic 
strip  is  expressed  in  the  classical  form. 

i,  t  de  (t) 

"W-llcrW^.Tf— *  (2-23) 

For  a  step  strain  history,  equation  (2.15), 

•<*>  -  3  eo  W*)  (2-24) 

and  the  viscoelastic  stress  intensity  factor  of  equation  (2.21)  can 
be  expressed  in  an  alternate  form 

"?w  VF  eo  W*>  (2-25) 

Either  equation  (2.25)  or  (2.21)  can  be  used  for  data  reduction. 


28 


Data  Reduction 

For  each  different  level  of  e  there  will  be  a  different  crack 

o 

velocity  history  produced.  Unfortunately  it  is  not  constant  in  time 
because  the  stress  in  equation  (2.21)  is  not  constant.  Parenthetically 
it  should  be  stated  here  that  if  a  means  for  experimentally  applying  a 
constant  stress  boundary  condition  could  be  employed,  then  the  results 
of  the  first  case  of  Appendix  B  could  be  used  and  nothing  would  vary 
with  time  in  the  strip.  For  a  given  level  of  stress  the  stress  intensity 
factor  would  remain  constant,  since  it  is  independent  of  crack  length; 
consequently  the  crack  velocity  would  be  constant,  and  it  would  be  a 
very  desirable  test.  In  fact  it  would  be  just  like  the  elastic  case 
for  which  this  test  has  been  used  by  Knauss  for  a  number  of  years 

Returning  to  data  reduction,  the  stress  and  crack  length  in  the 
strip  are  recorded  as  a  function  of  time,  as  shown  in  Figure  2-1+ . 


Figure  2-k     Representative  Data  From  Strip  Test 

From  the  crack  history,  the  crack  velocity  history  can  be  constructed  by 
taking  the  slopes  graphically  (Fig.  2-5).   A  cross-plot  can  then  be 
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made  of  stress,  or  equivalently  stress  intensity  factor,  versus  crack 
velocity.  This  is  the  master  curve,  or  fracture  characterization 
curve  shown  in  Figure  2-5.  Normally  these  tests  will  be  conducted  over 
a  range  of  temperatures  in  order  to  obtain  the  whole  curve. 


Figure  2-5 
Master  Fracture  Characterization  Curve 

I"l3l 
HertzlerL  J  has  obtained  such  characterization  data  over  a  range 

of  temperatures  from    15  F  to    85  F  for  a  CTFB  propellant.  He 
found  that  the  data  shifted  very  well  to  make  a  master  curve,  and  the 
shift  was  exactly  the  same  as  that  obtained  by  other  investigators  for 
the  relaxation  modulus.  His  results  are  shown  in  Figure  2-6,  Part 
of  the  scatter  in  the  data  is  attributed  to  the  fact  that  the  approximate 
expression  was  used  to  reduce  it. 

To  our  knowledge  this  is  the  first  viscoelastic  fracture  characteriza- 
tion of  propellant.  Some  scatter  in  the  data  is  caused  by  local  inhomo- 
geneitie6  in  the  material,  which  in  turn  cause  the  crack  velocity  to 
fluctuate  as  the  tip  proceeds  through,  or  around,  these  regions. 
Nevertheless  the  master  curve  is  well  defined  and  exhibits  no  more 
scatter  than  is  normally  noted  in  propellant  data  of  other  varieties, 
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even  modulus  data,  which  should  evince  the  smallest  amount  of  scatter 
of  all. 

The  master  curve  indicates  two  limit  values  that  are  important. 
The  low  velocity  region  gives  evidence  of  a  lower  bound  of  stress 
intensity  factor  for  which  there  will  be  no  crack  growth.  This  figure 
of  K..   =10  psi  -Tin  is  a  useful  guide  in  making  limit  checks  on  behavior 
of  cracks  in  rocket  motors,  from  which  their  structural  stability  can  be 
determined.  Any  loading  below  this  limit  should  be  safe,  at  least  for 
the  present  time;  however,  the  tendency  of  this  characterization  to  change 
with  time,  as  in  a  storage  condition,  is  not  known.  Future  testing 
could  be  directed  toward  this  end;  however  it  would  require  considerable 
amounts  of  propellant,  since  many  rather  large  specimens  are  required. 

Specimens  selected  for  this  characterization  were  marginally  small 
in  size  and  were  6"  x  1"  x  l/k"   and  bonded  to  wooden  tabs  along  the 
six  inch  dimension.    Five  different  specimens  were  ru»  for  each  of 
the  six  temperatures.  Of  course  it  would  be  best  if  the  entire  set  of 
tests  were  duplicated  at  each  time  interval  of  the  period  selected. 
From  this  information  a  good  description  of  fracture  behavior  over  the 
desired  shelf  life  could  be  obtained. 

In  lieu  of  that  data,  one  might  assume  on  the  basis  of  the  results 
in  Section  1  that  the  stress  intensity  factor  will  degrade  in  the  same 
manner  as  the  ultimate  strength.  This  quantity  is  routinely  measured  in 
biaxial  endurance  data  and  forms  the  basis  for  extrapolating  to  the 
storage  life  required.  With  this  information,  the  behavior  of  cracks 
existing  in  rocket  motors  that  are  currently  stored  could  be  predicted, 
and  the  reliability  of  the  missile  from  this  standpoint  could  be  assessed. 
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The  other  end  of  the  spectrum  is  also  of  interest.  It  appears 
that  if  the  data  were  extended  in  the  high  velocity  range,  the  ultimate 
limit  of  crack  velocity  would  be  approximately  10  inches  per  minute.  This 
is  indeed  slow  and  substantiates  the  earlier  assumptions  that  neglected 
the  inertial  terms  in  the  equilibrium  equations.  An  ultimate  limit 
strongly  influences  the  crack  history.  For  instance  consider  an  infinite 
sheet  with  a  central  crack  subjected  to  a  step  stress  at  infinity. 
From  equation  (2.10) 

K^t)  =  aoh(t)vficTtJ  (2.26) 

From  Figure  2-6,  if  K.(o)  is  greater  than  10  psi  \/in,  the  crack  will 
begin  to  grow.  As  it  grows  K..(t)  increases,  and  in  turn  the  crack 
propagates  faster.  This  process  continues  until  the  crack  reaches  a 
limit  velocity,  at  which  time  it  will  continue  to  propagate  across  the 
infinite  sheet  at  constant  speed,  even  though  the  stress  intensity  factor 
continues  to  increase. 

The  time  to  reach  the  constant  velocity  portion  of  the  history  is 
a  function  of  the  magnitude  of  a   .   If  a  were  large  enough,  the  crack 
could  begin  its  propagation  in  the  constant  velocity  region,  in  which 
case  the  behavior  would  be  elastic,  much  like  most  observed  behavior  with 
metals  or  glassy  polymers.  On  the  other  hand,  infinite  time  could  be 
required  if  K(0)  is  below  10  psi  Jiri   .  However  any  TC(0)  above  10  psi  J±ri 
eventually  will  produce  the  constant -maximum  velocity  propagation  for 
this  kind  of  geometry,  where  K  (t)  is  continually  increasing. 
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Solution  Method 

The  method  for  predicting  the  crack  behavior  in  viscoelastic 
materials  essentially  follows  the  method  proposed  by  Knauss     for 
Solithane^with  some  variations  for  the  relaxation  behavior.  For 
discussion  purposes,  consider  a  rocket  motor  with  a  bore  crack,  where 
it  is  desired  to  predict  what  the  crack  will  do  under  a  given  loading 
program. 

If  an  analytic  expression  for  the  stress  intensity  factor  of  the 
elastic  motor  geometry  can  be  obtained,  then  by  the  correspondence 
principle  previously  discussed,  the  viscoelastic  stress  intensity  factor 
is  obtainable.  Suppose  it  has  the  form 

K^(t)  =  \[c(t)]a(t)^TtT  (2.27) 

From  the  fracture  characterization  curve  similar  to  Figure  2-6, 

K^(t)  =  f(c)  (2.28) 

Inserting  equation  (2.28)  into  (2.27),  a  first  order  differential 
equation  in  crack  length  is  produced,  which  when  integrated  will  give 
the  desired  history.  Comparison  of  that  history  with  storage  life, 
or  burning  rates  under  operational  conditions,  will  determine  the 
detrimental  effect,  if  any,  that  the  crack  may  contribute. 

Often  times  motor  geometries  are  too  complex  to  allow  analytic 
expressions  for  the  stress  intensity  factor.   In  this  case  numerical 

methods  have  been  developed  to  find  elastic  stress  intensity  factors. 

("191 
Deverall  and  Lindsey     have  suggested  some  simple  methods  by  which 

these  can  be  obtained.  Subsequently  numerical  transforms  and  inversions 

would  be  necessary  to  obtain  the  final  solutions.  Because  this  is  an 
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involved  solution  process,  an  approximate  method  is  being  developed  and 
studied,  which  looks  very  promising.  It  is  based  upon  an  observation 
by  Dr.  Russell  Westman  of  UCLA,  that  the  local  stress  in  bodies  without 
cracks  can  often  be  related  to  the  stress  intensity  factor  for  the  same 
body  with  a  crack.  Progress  on  this  development  will  be  given  in  a 
later  report. 

Conclusions 

It  is  felt  that  the  viscoelastic  crack  behavior  for  materials  in 
which  the  crack  propagates  slowly  is  predictable  from  the  concepts 
developed  in  this  section.  Limitations  on  loading  and  geometry, 
however,  must  be  such  that  symmetric  crack  behavior  is  produced.  Further 
investigations  will  be  made  to  ascertain  if  the  concepts  developed  in 
Section  1  for  mixed  mode  propagation  can  be  applied  to  viscoelastic 
fracture.  Further  evaluations  of  the  validity  of  these  concepts  will 
be  made  using  existing  data  from  the  literature  and  in  generating  new 
experimental  crack  propagations  results  where  needed. 
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SECTION  3 
DEWETTING 


SECTION  3 

DEWETTING 
Introduction 

Propellant  differs  from  classical  viscoelastic  materials  through 
dewetting.   Since  dewetting  has  a  pronounced  influence  on  the  stress- strain 
law  for  strain  levels  experienced  in  many  rocket  motors,  it  is  necessary 
to  investigate  the  influence  of  dewetting  on  fracture. 

r  -I  £-   -I  £.       "I  rj~\ 

In  previous  reports  Lindsey   '       and  co-workers  have  developed 
an  isotropic  theory  for  dewettable  materials.  From  this  development  a 
means  has  been  formulated  by  which  stress  analyses  can  be  carried  forth 
while  accounting  for  dewetting.  In  that  theoretical  treatment  it  has  been 
postulated  that  a  stress  criterion  for  dewetting  exists  and  that  is  is 

expressible  in  terms  of  a  surface  in  principle  stress  space.  Furthermore 

Cl8] 
by  application  of  Drucker's     Fundamental  Theorem  for  stable  materials, 

[19] 
Wood     has  obtained  the  general  expression  that  describes  the  dewetting 

criterion  for  arbitrary  stress  states 


J2  +  a   I  =  Constant  (3.1) 


where 


9K2   dG 

In  order  to  evaluate  these  results  and  the  associated  assumptions, 
experimental  tests  were  developed  and  conducted  for  dewetting  in  uniaxial 
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and  triaxial  stress  fields. 

The  experiments  employed  a  liquid  dilatometer,  both  pressurized  and 
unpressurized,  and  the  poker chip  test  .   A  brief  synopsis  of  the  experi- 
ments is  presented  here.  At  a  later  date  the  detailed  treatment  of  the 
theory  and  experiments  will  be  presented  in  LCDR  J.  E.  Wood's  PhD  thesis. 
Dilatometer  Tests 

A  liquid- filled  dilatometer  using  Dow-Corning  DC-200  silicone  fluid 
was  employed  for  the  measurement  of  volume  change  in  uniaxial  tensile 
test  specimens.  Load  was  sensed  from  strain  gauges  mounted  on  a  flexure 
in  line  with  the  test  specimen  within  the  dilatometer  chamber.  Volume 
change  was  determined  by  visually  measuring  the  fluid  height  in  a  capillary 
tube  which  permitted  a  minimum  resolution  of  dilatometric  strain  of  0.01  per 
cent.  In  order  to  minimize  the  effect  of  thermally- induced  variations, 
fluid  height  readings  were  compared  with  readings  from  a  reference  chamber 
which  duplicated  the  test  chamber  in  all  thermal  characteristics.  Maximum 
test  duration  was  limited  to  thirty  minutes. 

A  number  of  different  formulations  of  highly  filled  composite  solid 
propellant  were  tested.  However  the  most  complete  set  of  tests  was 
carried  out  on  a  highly  filled  PBAN  propellant  supplied  by  United  Technology 
Corporation.  All  of  the  following  data  were  obtained  on  that  propellant. 
Constant  Load  Dilatometer  Tests 

Constant  load  uniaxial  tension  tests  were  conducted  in  the  dilatometer 
to  investigate  the  time- dependent  nature  of  the  onset  of  dewetting  under 
small  constant  loads.  Figure  3-1  shows  dilatation  vs.  time  for  a  series 
of  five  tests  of  the  UTC  PBAN.  The  rapid  increase  in  dilatation  observed  to 
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occur  soon  after  the  application  of  load  for  the  three  higher  stresses 
gives  a  clear  indication  of  dewetting.  Because  of  the  very  small 
amounts  of  volume  change  produced  in  these  tests  at  the  lower  stress 
levels,  it  is  difficult  to  determine  the  presence  or  absence  of  dewetting 
from  the  figure.  Further  amplifying  information  is  obtained  from  the 
variation  in  Poisson's  ratio  with  time  for  the  same  specimens,  shown  in 
Fig.  3-2.  In  this  plot  the  rapid  decrease  in  v   associated  with  dewetting 
is  seen  to  occur  immediately  at  55-^  psi  and  31-5  psi.  The  response  for 
the  three  lower  stresses  shows  that  the  drop  in  u  is  delayed  in  time  and 
reduced  in  magnitude  with  decreasing  stress. 

Qualitatively  similar  behavior  was  observed  for  all  composite  pro- 
pellant  tested,  including  PBAN,  PBAA,  CTPB,  and  HTPB  formulations. 

From  these  tests,  it  was  concluded  that  although  the  onset  of 
dewetting  in  the  propellants  tested  does  demonstrate  some  time  dependence 
for  sufficiently  small  values  of  stress,  the  small  amount  of  dilatation  which 
results  makes  evaluation  of  this  dependence  impractical  in  the  dilatometer. 
However,  a  reasonably  narrow  zone  of  stress  values  exists  above  which  dewetting 
effects  are  negligible,  at  least  for  loads  of  short  duration.  It  is  in  this 
sense  that  a  stress  criterion  for  dewetting  will  be  interpreted. 
Constant  Strain  Rate  Dilatometer  Tests. 

These  tests  were  carried  out  to  identify  the  stress  criterion  for  the 
onset  of  dewetting  in  uniaxial  tension  and  uniaxial  tension  with  superimposed 
hydrostatic  pressure  at  a  strain  rate  of  one  per  cent  per  minute,  for  which 
the  material  response  prior  to  dewetting  was  quasi- elastic. 

Figure  3-3  shows  the  stress- strain-dilatation  response  for  four  dupli- 
cate tests  of  PBAN  specimens.  The  dilatation  increment  plotted  has  the 
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initial  linear  response  subtracted  out,  thus  the  intercept  of  the  dilata- 
tion curve  with  the  axis  indicates  the  onset  of  dewetting. 

Following  the  above  tests,  pressurized  dilatation  tests  were  under- 
taken in  an  attempt  to  determine  the  onset  of  dewetting  under  a  combined 
state  of  stress.  Tests  were  run  with  superimposed  pressures  of  10,  20, 
and  U0  psig.  Maximum  test  pressure  was  limited  to  kO   psig  because  of 
excessive  creep  of  the  plexiglas  walls  of  the  dilataneter  chamber  experienced 
at  higher  pressures.  More  data  scatter  was  encountered  with  these  pressurized 
tests;  however,  no  significant  differences  from  the  unpressurized  tests 
either  in  the  stress  at  onset  of  dewetting  or  in  the  character  of  the 
response  could  be  detected. 
Pokerchip  Triaxial  Tension  Test 

In  order  to  investigate  dewetting  behavior  in  a  stress  state  approxi- 
mating hydrostatic  tension,  a  pokerchip  apparatus  was  employed.  For 
materials  with  Poisson's  ratio  near  0.5,  the  pokerchip  triaxial  tension 
test  provides  a  state  of  approximately  hydrostatic  tension  within  the 
central  region  of  the  test  specimen.  Using  the  dewetting  criterion  previously 
discussed,  it  is  found  from  the  analytical  solution  for  stresses  within  the 
pokerchip  that  the  dewetting  criterion  is  first  satisfied  at  the  center  of 
the  specimen,  within  the  region  of  hydrostatic  tension.   Thus  the  poker- 
chip test  provides  a  means  of  initiating  dewetting  under  a  known  triaxial 
tension  state. 

In  order  to  interpret  the  results  of  the  pokerchip  test  properly, 
it  was  necessary  to  determine  the  influence  of  dewetting  on  experimentally 
measurable  quantities.  To  investigate  this  question  a  computerized  finite 
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element  analysis  of  the  pokerchip  was  performed.  The  Rohm  and  Haas 
AMG  032A  structural  analysis  program  for  axisymmetric  structures  and 
loads  was  adapted  to  the  IBM  36O-67  computer  at  the  Naval  Postgraduate 
School.  This  program  is  designed  to  handle  materials  with  Poisson's 
ratio  at  or  near  0.5,  and  allows  variation  of  material  properties  from 
element  to  element,  thus  it  was  suitable  for  the  task. 

To  test  the  program  and  establish  convergence,  the  program  was 
used  to  analyze  a  pokerchip  specimen  loaded  with  a  uniform  axial 
elongation.  Results  were  compared  with  the  approximate  analytical 
solution  given  by  the  LSWZ  c       analysis.  The  geometry  chosen  was  that 
which  would  be  used  in  later  testing  -  a  disk  V  in  diameter  and  0.25" 
thick.  Various  modulus  values  representative  of  composite  solid  pro- 
pellant  were  investigated,  with  Poisson's  ratios  in  the  range  0.^9  to 
0.^99*  The  program  was  found  to  perform  satisfactorily,  giving  close 
correlation  with  the  LSWZ  solution.  The  program  was  designed  to  have  the 
capability  of  handling  incompressible  material  (v  -   .5),  however,  it  was 
found  that  for  the  very  stringent  displacement  boundary  conditions  imposed 
by  the  pokerchip  test,  convergence  became  less  easy  to  achieve  as  u 
approached  0.5.  By  utilizing  the  maximum  allowable  number  of  elements  in 
the  radial  direction,  it  was  possible  to  obtain  good  solutions  for  x>  -   0.^99 
Table  3-1  gives  a  comparison  of  stresses  calculated  by  the  program  and  by 
the  LSWZ  solution  for  this  value.  The  comparison  is  quite  good,  diverging 
slightly  only  near  the  outer  radius  where  the  LSWZ  solution  posseses 
singularities. 
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To  model  the  effect  of  dewetting  in  the  pokerchip  test,  it  was 
considered  that  the  initial  effect  of  dewetting  would  be  the  alteration 
of  material  properties  in  a  small  region  at  the  center  of  the  specimen. 
For  the  finite  element  representation,  a  two-phase  specimen  was  employed, 
having  an  inner  core  l/2"  in  diameter  to  which  reduced  values  of  E  and  v 
were  assigned,  while  the  properties  in  the  rest  of  the  pokerchip  were 
unmodified.  Figure  3-h   shows  a  cross-section  of  the  finite  element  grid 
employed. 
Results 

It  was  found  that  most  aspects  of  the  solution  for  the  two-phase 

specimen  showed  only  slight  changes  as  properties  within  the  core  were 

altered.  But  a  significant  change  was  noted  in  the  state  of  stress 

around  the  center  of  the  specimen,  especially  within  the  core  itself, 

for  slight  variations  in  Poisson's  ratio  in  the  core  material.  Figure 

3-5  shows  the  computed  axial  stress  at  the  midplane  versus  radius  for 

various  values  of  v   at  the  central  core.  For  a  reduction  of  u    of 

core 

0.8$  (from  .1*99  "to  «^95)  the  stress  at  the  center  drops  by  approximately 
10$,  with  progressively  greater  reductions  occurring  as  u     is  further 
decreased.  The  indicated  percentage  change  in  the  stress  state  at  the 
center  is  at  least  an  order  of  magnitude  greater  than  that  occurring  in 
any  of  the  customarily  measured  quantities  such  as  total  load  or  radial 
displacement  on  the  periphery. 

On  the  basis  of  this  analysis  it  appears  that  a  pressure  transducer 
located  at  the  center  of  one  of  the  platens  of  the  test  rig  could 
easily  detect  a  deviation  in  the  local  stress  state  due  to  a  small 
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TABLE  3-1     Midplane   Stresses  in  Pokerchip  Test  Specimen 


Specimen  Diameter:  1+.0",  thickness  0.25" 

Material  Properties:     E  =  300  psi,   u  =  0.1+99 
Loading:     Uniform  axial  elongation  -  Az/t  =  0.005 


Radius 

Finite  Element  Sol. 

LSWZ   Sol. 

Finite  Element  Sol. 

LSWZ   Sol. 

Inches 

a 
r 

CTe 

r          9 

a 
z 

a 
z 

PSI 

PSI 

PSI 

PSI 

PSI 

.0625 

7^.5 

7^.5 

7I+.I 

76.2 

75.1+ 

.1875 

73.8 

73.8 

73.6 

75.7 

75.0 

.3125 

72.8 

72.8 

72.5 

7^.7 

73.9 

.*+375 

71.1 

71.2 

70.9 

73.0 

72.3 

.5625 

69.I 

69.I 

68.8 

71.0 

70.1 

.6875 

66.3 

66.3 

66.0 

68.2 

67.5 

.8125 

63.1 

63.1 

62.9 

65.O 

61+. 2 

•  9375 

59.1 

59-1 

59-0 

61.0 

60.1+ 

1.0625 

5U.8 

5^.8 

5U.6 

56.7 

56.0 

1.1875 

1+9-5 

1+9.6 

49.6 

51.5 

51.0 

1.3125 

1+3.9 

¥+.0 

1+3.9 

45.9 

1+5.3 

1.1*375 

37.  k 

37.5 

37.6 

39.5 

39.0 

1.5625 

30.3 

30.1+ 

30.5 

32.1+ 

32.0 

1.6875 

22.5 

22.6 

22.8 

21+.5 

21+.2 

1.8125 

13.7 

1I+.2 

14.3 

16.6 

15.8 

1.9375 

3.8 

3.5 

5.0 

5.1 

6.5 

1+6 


alteration  of  properties.  While  at  present  the  extent  to  which 
dewetting  in  hydrostatic  tension  does  alter  the  Poisson's  ratio 
has  not  been  determined,  such  an  installation  would  appear  to  be  at 
least  as  sensitive  to  Poisson's  ratio  as  the  dilatometer  is  for  uni- 
axial testing,  and  should  provide  an  equally  good  indication  of  dewetting 
initiation. 

Another  finite  element  problem  was  worked,  which  provided  for  a 
section  of  flexible  boundary,  just  above  the  dewetted  core,  where  the 
pressure  transducer  would  be  mounted  in  the  otherwise  rigid  grip.  This 
analysis  was  conducted  to  determine  if  introduction  of  the  pressure 
transducer  affected  the  stress  field  in  the  specimen.  Even  by  allowing 
the  displacement  at  the  boundary  to  be  200  times  the  expected  displacement 
of  the  transducer  diaphragm,  the  change  in  the  stress  field  was  imperceptable , 

Accordingly  a  pressure  transducer  with  a  l/2"  diameter  contact  area 

was  imbedded  in  the  center  of  one  of  the  k"   diameter  pokerchip  platens, 

[21] 
similar  in  construction  to  that  reported  by  Harbert    .  Figure  3-6 

shows  the  stress- strain  plot  obtained  from  one  of  the  pokerchip  tests  run 

at  an  axial  strain  rate  of  0.007/min.  The  dashed  lines  indicate  the  elastic 

solution  for  the  specimen  with  E  =  520  psi  and  u  =  .U-95  •     At  an  axial  strain 

of  0.0075 i   the  stress  sensed  by  the  pressure  transducer  at  the  center  of  the 

specimen  shows  a  sharp  reduction  from  the  initial  linear  slope.  The  stress 

level  at  this  point,  9^  psi  in  this  case,  is  taken  as  the  point  of  dewetting 

initiation. 

Pokerchip  tests  instrumented  with  a  central  pressure  transducer  have 

been  reported  by  Harbert.  For  his  tests  performed  at  a  constant  rate  of 
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axial  elongation,  a  sharp  break  in  the  pressure  transducer  reading  was 
observed  at  a  characteristic  level  of  stress  which  was  identified  as 
the  point  of  initial  failure.  From  the  analysis  reported  above,  it 
appears  that  such  a  result  may  indicate  only  a  slight  reduction  in 
load- carrying  capability  which  may  be  at  a  stress  well  below  the  actual 
failure  stress. 
Comparison  of  Experimental  Data  with  Theoretical  Criterion 

The  theoretically  derived  dewetting  criterion  referred  to  previously 

is 

2 
J„  +  al   =  Constant  (3.1) 

in  which  J0  =  is. .S. .  is  the  second  invariant  of  the  deviatoric  stress 

tensor  and  I-|=avi,  is  the  first  invariant  of  the  stress  tensor,  while 

2 
G    dK 

°  =  9^      dG 

is  a  material  parameter  dependent  upon  the  change  in  the  shear  and 
bulk  moduli  of  the  material  as  dewetting  initiates.     From  data  for  the 
unpressurized  uniaxial  tensile  test  runs  shown  in  Fig.   3-3  ,   ot  was  calcu- 
lated to  be  0.0026,  making  the  dewetting  criterion  for  this  material  at 
the  temperature  and  rate  tested: 

J2  +  0.0026  1^    =     180 

2 
For  the  uniaxial  tensile  tests   j_  =  a  /3   and   I..  =  o    -  3P  . 

2    u'  1    u 

For  the  pokerchip  triaxial  tests  J?  =  0     and   I,  =  3CT+  . 

The  stresses  measured  at  the  point  of  dewetting  in  the  various  tests 
described  are  plotted  in  Fig.  3-7.  The  solid  curve  shows  the  theoretical 
dewetting  criterion. 
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While  the  limited  amount  and  range  of  the  data  presented  cannot 
alone  confirm  the  validity  of  the  stress  criterion,  it  is  satisfying 
to  note  the  good  correlation  with  the  triaxial  tension  results. 
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SECTION  k 
CONCLUSION 


SECTION  k 
CONCLUSION 

Several  aspects  of  polymer  fracture  have  "been  discussed  in  this 
report.  A  theory  for  crack  instability  under  mixed  mode  conditions  has 
"been  proposed,  and  to  a  limited  extent  verified,  for  a  glassy  polymer. 
Normally  advanced  concepts  such  as  this  are  first  attempted  on  the 
simpler  elastic  materials  where  there  exists  a  higher  probability  of 
success.  Then  once  a  foothold  has  been  established,  it  can  be  generalized 
and  made  applicable  to  viscoelastic  materials.  The  developments  and 
results  in  Section  1  constitute  work  of  that  type  and  further  work  will 
be  directed  toward  generalizing  it  to  he  applicable  to  viscoelastic 
materials  and  propellants. 

Viscoelastic  materials,  being  more  complex,  do  not  have  so  advanced 
a  fracture  theory  as  do  elastic  materials  *  The  evolution   of  theory  lags 
that  of  elastic  materials  but  is  following  along  in  a  very  similar  pattern 
of  development.  The  concept  of  stress  intensity  factor  for  viscoelastic 
materials  is  used  very  extensively  in  this  development,  much  more  so  than 
in  the  current  literature  on  the  subject.  The  methods  proposed  by  which 
materials  may  be  characterized  for  fracture  and  later  analyzed  for 
fracture  are  still  hypothetical  in  nature  and  must  be  evaluated 
experimentally.  This  work  is  presently  being  pursued. 

An  alternate  method  of  development  based  upon  an  energy  approach  as 
performed  by  Knauss  and  his  coworkers,  leads  by  a  very  complex  analysis 
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to  a  dependence  in  the  solution  on  an  undetermined  length  parameter, 
which  must  be  evaluated  experimentally.  The  approach  outlined  in 
Section  2  avoids  this  problem  by  going  directly  to  the  stress  intensity 
factor,  and  it  becomes  much  more  simple  mathematically. 

Finally  propellant  fracture  to  be  complete  must  incorporate  dewetting 
effects.  A  theory  to  do  this  in  quasielastic  situations  has  emerged  and 
has  been  experimental  1 y  evaluated  for  a  limited  number  of  stress  states. 
Considerable  work  has  been  expended  in  developing  this  theory  and  LCDR 
J.  E.  Wood  has  completed  a  PhD  research  project  in  this  area.  His 
thesis  is  now  in  preparation  and  will  treat  this  topic  in  detail.  Copies 
will  be  forwarded  when  available.  A  brief  synopsis  of  just  a  few  of  the 
re  suits  are  included  in  Section  3  to  provide  an  indication  of  the  type  of 
work  that  has  been  done. 

An  evaluation  of  what  analyses  can  be  performed  at  this  time  leads 
to  the  following  conclusions: 

1.  Predictions  of  crack  instability  for  purely  elastic  materials 
can  be  made  for  two  dimensional,  through  the  thickness  cracks 
under  arbitrary  loading. 

2.  Predictions  of  crack  history  can  be  made  for  purely  viscoelastic 
materials  with  no  dewetting  for  two  dimensional,  symmetric 
geometries  and  loads. 

3.  Dewetting  can  be  incorporated  into  the  stress  analysis  for 
quasielastic  loading  conditions  such  as  cool -down.  Study  of  this 
problem  is  presently  underway  and  the  associate  problem  of  the 
temperature  dependence  of  the  dewetting  criterion  is  also  being 
studied.  Measurements  of  dilatation  over  a  range  of  temperatures 
has  been  completed.  The  step  to  incorporate  dewetting  into 
viscoelastic  materials  has  not  yet  been  made. 
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APPENDIX  A 

EXTENDED  CORRESPONDENCE  PRINCIPLE 

[a] 
The  extension  of  the  correspondence  principle  by  Graham    to 

cover  crack  propagation  problems  is  a  key  step  in  the  development  of 

viscoelastic  fracture  as  it  is  done  in  this  report.  As  a  result,  a 

synopsis  of  the  central  theme  of  the  paper  is  included  here;  however 

the  development  does  not  parallel  that  of  the  paper  but  enlarges  upon 

it.  It  is  hoped  that  by  reading  both  this  appendix  and  the  paper,  a 

clear  view  of  the  basis  for  Section  2  of  this  report  will  be  obtained. 

Definition  of  the  problem 

The  correspondence  principle  is  clearly  applicable  whenever  the 

type  of  boundary  condition  prescribed  is  the  same  at  all  points  of  the 

boundary.  For  mixed  boundary  value  problems  (i.e.,  problems  for  which 

different  field  quantities  are  prescribed  over  separate  parts  of  the 

boundary)  the  method  is  still  applicable  provided  the  regions  over 

which  different  types  of  boundary  conditions  are  given  do  not  vary 

with  time.  (We  are,  of  course,  assuming  that  the  region  occupied  by  the 

body  does  not  vary  with  time.)  There  remain  those  viscoelastic  mixed 

boundary  value  problems  where  the  regions,  over  which  different  types 

of  boundary  conditions  are  given,  do  vary  with  time.  Particular 

examples  are  indentation  and  crack  propagation  problems.  For  problems  of 

this  type  there  will  be  points  of  the  boundary  at  which  only  partial 

histories  of  some  field  quantities  will  be  prescribed.  When  this  is  the 

case  the  transforms  of  these  quantities  are  not  directly  obtainable  and 

the  classical  correspondence  principle  is  not  applicable. 
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We  desire  the  solution  to  the  viscoelastic  problem  shown 
diagramatically  in  Figure  A-l. 


Crock    Surface 


Ltr\e   o"f  propagation 


FIGURE  A-l 
Viscoelastic  Crack  Problem 


In  figure  A-l,  the  material  below  the  line  of  crack  propagation 
has  been  removed  and  subsequently  replaced  with  the  appropriate 
boundary  conditions,   a  (x,t),b(x,t)   andc(x,t)  represent  the  elements 
of  any  column  of  the  matrix     in  A.l. 


ooaauuuu 
ssnnssnn 


uouaauau 
nnssnnss 


auouuouo 
nnssnnss 


(A.l) 


u  and  u  denote  normal  and  tangential  components  of  the  boundary 
displacements,  and  o     and  o     represent  similar  quantities  for  stress. 
The  boundary  conditions  for  the  problem  are: 

aV(x, t)  =  AV(x,t)   on   B 

bv(x,t)  =  Bv(x,t)   on   Bx(t)  (A. 2) 


cV(x,t)  =  0 


on 


B0(t) 
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where  B  (t)  +  B2(t)  =  B,    B,(t)  is  monotonically  increasing,  and  B 
is  the  boundary  of  a  fixed  region  R. 

The  solution  is  obtained  by  first  solving  an  elastic  problem 
with  changing  boundaries  as  shown  in  Figure  A-2. 


e        e        e 

a.  .(x,t),  e.  .(x,t),  u.(x,t)  are  solutions  in  R. 


Q*M 


B^ti 


be5? 


5r(t)  


FIGURE  A-2 
Equivalent  Elastic  Problem 

g 

With  the  solution  to  this  problem,  the  boundary  values  of  b  (x,t) 

on  Bp(t)  are  computed,  and  the  values  of  c  (x,t)  are  computed  for  B  (t) 

With  these  values,  now  define 


be(x,t)  =  Be(x,t) 


on   B 


(A. 3) 


ce(x,t)  =  k(E,u)  Ce(x,t)   on   B 


The  constraints  on  the  problem  must  be  such  that  the  new  boundary 
conditions  B  (x,t)  are  independent  of  material  properties.  This  is 
necessary  for  the  classical  correspondence  principle  to  work  in  the  next 
step  of  derivation. 

We  can  now  define  the  same  elastic  problem  in  terms  of  new 
boundary  conditions  as  shown  in  Figure  A-3 
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e        e        e 

tf..(x,t),  e  (x,t),  ut(x,t) 


ari^t)      b'OM:) 


FIGURE  A- 3 
Elastic  Problem  with  New  Boundary  Conditions 
Because  of  the  nature  of  the  new  boundary  conditions,  transforms  can  be 
computed  for  them  since  they  are  now  defined  for  all  time.   That 
condition  is  shown  in  Figure  A-k. 


*    *    * 
o\  .,  e.  .,  u.   are  solutions  in  R 


B.a) 


CL«*   b«* 


e^ce 


FIGURE  A- k 


Corresponding  Elastic  Problem 


This  is  the  corresponding  solution  for  the  viscoelastic  problem 
shown  in  Figure  A- 5. 
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CT..(x,t),  e..(x,t),  u.(x,t)  are  viscoelastic  solutions  in  R 


ofc,0 


bMJ 


■&#)■ 


av6t,t)        bvfr>t) 


V*> 


FIGUKE  A-5 
Viscoelastic  Problem  Equivalent  to  Figure  A-l 


Therefore  the  correspondence  principle  is  given  by  the 
classical  relationships. 


v*  r  e*     .  ;>-, 

*ii(x^  =    LGi.i(x>p)  J 


e* 

.( 


E  -  E* 


v*  e*     -, 

:..(x,p)  =   Leid(x>p)  J 


E  -  E* 


(A.U) 


v*         p  e*      -, 
u.  (x,p)  =   [u._  (x,p)  J 


E  -»  E* 


This  extended  version  of  the  correspondence  principle  is  used 
to  generate  theoretical  developments  pertaining  to  viscoelastic  fracture 
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Extension  of  Graham's  Paper 

An  extension  of  Graham's  development  can  be  made  by  introducing 
a  third  boundary,  B_,  which  is  constant  in  time;  so  that  now  the 
unchanging  body  region,  R,  is  bounded  by 

B  =  B1(t)  +  B2(t)  +  B3 
Since  B  is  constant, 

B.,(t)  +  B  (t)  =  constant, 
as  in  Graham's  paper. 

B_  can  consist  of  an  arbitrary  number  of  segments,  upon  each  of  which 
the  boundary  conditions  may  be  different;  however  upon  none  of 
the  segments  can  the  type  of  boundary  condition  change  with  time. 

Such  a  difinition  of  the  boundary  conveniently  facilitates  the 
application  of  the  correspondence  principle  to  finite  geometries. 
For  instance  consider  the  infinite  strip  geometry  with  displacement 
boundary  conditions  as  shown  in  Figure  A-6. 


B= 


I  l  I  I  i   /  l   I   I    //////  / 


E>z(t) 


Figure  A-6  Infinite  Strip  Subjected  to  Displacement  Boundary  Conditions 
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The  boundary  conditions  on  B_  are  known  for  all  time,  and  therefore 
their  transforms  can  be  computed.  Now  applying  the  original  methods 
of  Graham,  the  entire  development  goes  through  as  before;  however  now 
a  more  versatile  geometry  is  available  for  which  the  correspondence 
principle  is  applicable.  The  prime  limitation  upon  generating  solutions 
is  the  requirement  that  B  (x,t)  in  equation  (A. 3)  be  independent  of 
material  properties. 
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APPENDIX  B 
STRESS  ANALYSIS  OF  AN  INFINITE  STRIP 

Viscoelastic  crack  propagation  has  been  extensively  studied  by 

["10  12.     22  23  2U~1 
Knauss  and  coworkers.   '   '   '        '      Much  of  the  work  has  been 

done  on  an  infinite  strip  with  a  semi  infinite  crack.  A  simplified, 

but  oft  times  adequate,  solution  method  is  given  in  this  appendix  for 

the  stress  intensity  factor  of  the  strip.  The  elasticity  theory 

required  is  very  elementary,  and  it  is  presented  here  for  completeness 

in  the  development  of  the  viscoelastic  stress  intensity  factor. 

Elastic  Stress  Analysis  of  Strip 


I    t    t    I    t   I   t   t    II   It   t   t 


I    I    I    I    I    I    I    I   I   I   I    I  I    j 


Figure  B-l  Stress  Boundary  Conditions 
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T    =  0 


The  Stress  Boundary  Conditions  are  specified  as 

y  =  ±b  °y  =  CT 

No  x  =  constant  boundaries  (infinite) 
Zero  stress  in  the  x  direction  at  infinity. 
Consider  a  stress  function  of  the  form 

0  =  -  x 
Obviously  this  function  satisfies  the  biharmonic  equation.  The 
components  of  stress  produced  by  this  function  are: 

=  0 


(B.l) 


a  =  0 
x 


ay  -a 


xy 


(B.2) 


This  stress  field  satisfies  the  boundary  conditions  and  therefore 
constitutes  the  solution  to  the  problem.  The  corresponding  strain  field 
is: 

If  a  semi -infinite  crack  is  inserted  in  the  strip,  as  shown  in 
Figure  B-2,  the  solution  obtained  above  will  be  valid  out  in  front  of 
the  crack  at  some  distance  where  the  disturbance  caused  by  the  crack 
has  decayed  to  zero. 


.t   I    I    I   I   M  M   1  I   II 


'   •   M    I    I    I    J   I    II   I   I   11 


Figure  B-2  Infinite  Strip  with  a  Seminfinite  Crack 
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To  the  far  left,  the  sheet  is  completely  unloaded  and  stress  free.  To 

the  far  right  it  is  under  constant  stress  as  if  there  were  no  crack, 

and  in  the  center  there  is  a  transition  region  of  unknown  character 

between  the  two  constant  fields.  If  the  crack  advances  by  an  amount  Ac, 

the  transition  region  moves  ahead  Ac.  The  constant  stress  field  is 

reduced  in  dimension  by  £c,  and  the  unstressed  region  is  increased  by 

the  same  amount.  The  strain  energy  lost  from  the  constant  stress  field 

by  the  crack  moving  ahead  Ac  is 

AU  =  ^  a   e  2btAc  (B.k) 

2  y  y  v   ' 

where  t  is  the  thickness  of  the  sheet.   In  the  limit  of  infinitesimal 
increments  in  crack  length, 

%   =  a  e  bt  =  g_bt  (B.5) 

dC     y  y      E  v    ' 

From  Irwins  development,  which  related  stress  intensity  factor  and 

strain  energy  release  rate  for  cracks  in  plane  stress, 

Kl  -  1   dU  (n  ^ 

T'tHS  (B-6) 

Thus  the  stress  intensity  factor  for  the  infinite  strip  is 

Kx  =  Ov/b  (B.7) 

It  is  of  interest  to  note  that  this  geometry  is  independent  of 

crack  length  and  therefore  possesses  a  constant  stress  intensity  factor 

as  the  crack  propagates  across  the  specimen.  For  this  reason  this 

geometry  lends  itself  well  to  laboratory  measurements.  Mueller1 

has  obtained  a  solution  of  a  pressurized  crack  in  a  semi- infinite  strip, 

and  from  his  results,  one  can  determine  the  length  of  crack  required  to 

produce  the  same  effect  as  an  infinite  strip. 
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The  difficulties  encountered  in  experimentally  holding  a  constant 
stress  at  the  boundary  as  the  crack  moves  is  great  enough  to  pursuade 
most  investigators  to  use  displacement  boundary  conditions. 

Displacement  Boundary  Conditions 

For  the  same  strip  geometry  subjected  to  displacement  boundary 

conditions  on 

y=±b  u=0         v=v  (B.8) 

o 

Guess  a  displacement  field  that  is  everywhere  zero  for  u  and  is  odd 

in  y  for  v. 

u(x,y)  =  0  (B.9a) 

v(y)  =  aQ  +  axy  +  a^3  ...  (B.$b) 

On  the  center line  of  the  sheet,  v  must  be  zero  by  symmetry;  therefore 

a  =0.  The  strain  field  corresponding  to  these  displacements  is 

2 
e  =  0       e  =  a..  +  3a_y  +  ...   e  =e  (x,v)     Y   =0 
x  y    1     3  z  zv  'J;      xy 

(B.10) 

The  stresses  that  result  from  such  boundary  conditions  are  calculated 

from  the  stress -strain  law, 

0  =  a     -   va  (B.lla) 

x     y 

e^  +   3a3y2  +  ...  =  |[a  -  vax]  (B.llb) 

e  =  ^a  +  a  ]  (B.llc) 

z    E  x    y 

From  (B.lla) 

o     =  vo   ,  (B.12a) 

x     y' 

and 

2 

^  +  3a3y2  +  . . .  =  [^>y  (B.lZb) 
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ez  =  -g(l  +  v)ay  (B.12c) 

The  two  dimensional  equilibrium  equations  simplify  to 

.da  da 

v  ^~  =  °  ^  =  0  (B.13) 

Using  equation  (B.12b),  equilibrium  is  satisfied  if 

E 
^[6a  y  +  .. .]  =0 

1-v 
Setting  a~  and  all  coefficients  of  higher  order  terms  to  zero,  the 
solution  reduces  to 

u(x,y)  =  0  (B.lUa) 

v(y)  =  s^y  (B.l4b) 

v 
The  boundary  conditions  are  satisfied  if  a.  =  — ;  therefore, 

v 
U(x,y)  =  0  v(y)  =  ^°-y  (B.15) 

The  strains  are  obtained  by  differentiation. 

v  v 

e=0       e   =  C2     «  ■  -v(l-v)^2    Y   =0        (B.16) 
x  y   b       z     v   'b     Txy  v    ' 

and  the  stress  field  is  expressed  by 


o 
x 


=  r_vEll£  a     -fJLJ^       a     =0  t        =0  (B.17) 

L,     2J  b  y      L,     2J  b  z  xy  " 

1-v  1-v 


As  before,  when  a  crack  is  inserted  into  this  specimen  the  constant 
stress  field  of  equation  (B.17)  is  present  at  a  distance  out  in  front 
of  the  crack.  As  the  crack  advances,  the  domain  of  constant  stress 
far  in  front  of  the  crack  is  reduced  by  Ac,  and  the  domain  of  zero 
stress  far  to  the  rear  of  the  crack  is  increased  by  Ac.  The  energy  lost 
in  the  process  of  the  crack  moving  by  Ac  reduces  to  the  same  expression 
as  before,  even  though  the  stress  and  strain  fields  are  different. 
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AU  =  £  a   e  2bt  Ac 
2  y  y 


In  the  limit  of  infinitesimal  crack  increments 

,v  v2 


dU    E 


dC   ,   2 
1-v 


(r) 


bt 


Once  again  from  Irwins  relationship  in  equation  (B.6) 


(B.18) 


(B.19) 


■'-jfetrt* 


or  from  (B.17) 


K 


l   y 


(B.20a) 


(B.20b) 
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